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Abstract. Solutions of a smooth first order dynamic equation can be made Lyapunov stable at 
will by the choice of an appropriate time-dependent Riemannian metric. 

1 Introduction 

The notion of the Lyapunov stability of a dynamic equation on a smooth manifold implies 
that this manifold is equipped with a Riemannian metric. At the same time, no preferable 
Riemannian metric is associated to a first order dynamic equation. Here, we aim to study 
the Lyapunov stability of first order dynamic equations in non-autonomous mechanics with 
respect to different (time-dependent) Riemannian metrics. 

Let us recall that a solution s(t), t G R, of a first order dynamic equation is said to 
be Lyapunov stable (in the positive direction) if for to G K and any e > 0, there is 5 > 
such that, if s'(t) is another solution and the distance between the points s(t ) and s'(t ) is 
inferior to 5, then the distance between the points s(t) and s'(t) for all t > to is inferior to e. 
In order to formulate a criterion of the Lyapunov stability with respect to a time-dependent 
Riemannian metric, we introduce the notion of a covariant Lyapunov tensor as generaliza- 
tion of the well-known Lyapunov matrix. The latter is denned as the coefficient matrix of 
the variation equation |], |||], and fails to be a tensor under coordinate transformations, un- 
less they are linear and time-independent. On the contrary, the covariant Lyapunov tensor 
is a true tensor field, but it essentially depends on the choice of a Riemannian metric. We 
show the following (see Propositions 0, and || below). 

(i) If the covariant Lyapunov tensor is negative definite in a tubular neighbourhood of 
a solution s at points t >t , this solution is Lyapunov stable. 

(ii) For any first order dynamic equation, there exists a (time-dependent) Riemannian 
metric such that every solution of this equation is Lyapunov stable. 

(iii) Moreover, the Lyapunov exponent of any solution of a first order dynamic equation 
can be made equal to any real number with respect to the appropriate (time-dependent) 
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Riemannian metric. It follows that chaos in dynamical systems described by smooth (C°°) 
first order dynamic equations can be characterized in full by time-dependent Riemannian 
metrics. 



2 Geometry of first order dynamic equations 

Let R be the time axis provided with the Cartesian coordinate t and transition functions 
t' = t+const. In geometric terms a (smooth) first order dynamic equation in non- 
autonomous mechanics is defined as a vector field 7 on a smooth fibre bundle 

tt:Y^R (1) 
which obeys the condition j\dt = 1, i.e., 

1 = d t + 1 k d k . (2) 
The associated first order dynamic equation takes the form 

t=l, V k = 1 k {t,yi)d k , (3) 

where (t,y k ,t,y k ) are holonomic coordinates on TY. Its solutions are trajectories of the 
vector field 7 (§). They assemble into a (regular) foliation T of Y. Equivalently, 7 fl2|) is 
defined clS db connection on the fibre bundle ([I]). 

A fibre bundle Y ([!]) is trivial, but it admits different trivializations 

Y = R x M, (4) 

distinguished by fibrations Y — » M. For instance, if there is a trivialization (^) such that, 
with respect to the associated coordinates, the components j k of the connection 7 @ are 
independent of t, one says that 7 is a conservative first order dynamic equation on M. 

Hereafter, the vector field 7 (0) is assumed to be complete, i.e., there is a unique global 
solution of the dynamic equation 7 through each point of Y. For instance, if fibres of 
Y — ► R are compact, any vector field 7 @ on Y is complete. 

Proposition 1. If the vector field 7 is complete, there exists a trivialization (fjj ofY 
such that any solution s 0/7 reads 

s a (t) = const., t e R, 

with respect to associated bundle coordinates (t,y a ). 
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Proof. If 7 is complete, the foliation T of its trajectories is a fibration ( of Y along these 
trajectories onto any fibre of Y, e.g., Y t=0 = M. This fibration yields a desired trivialization 
I- □ 

One can think of the coordinates (t, y a ) in Proposition [TJas being the initial date coordi- 
nates because all points of the same trajectory differ from each other only in the temporal 
coordinate. 

Let us consider the canonical lift Vj of the vector field 7 (p[) onto the vertical tangent 
bundle VY of Y — > R. With respect to the holonomic bundle coordinates (t,y k ,y k ) on 
VY, it reads 

V 1 = 1 + d jl k y j d k , d k = -^. 

This vector field obeys the condition V^f\dt = 1, and defines the first order dynamic 
equation 

t = l, y k = 1 k (t,f), (5a) 

V* = d j7 k (t,y l )yi (5b) 



on VY. The equation ( [5a]) coincides with the initial one (|3]). The equation fl5H5|) is the 
well-known variation equation. Substituting a solution s of the initial dynamic equation 
fl5a|) into (pbl), one obtains a linear dynamic equation whose solutions s are Jacobi fields of 



the solution s. In particular, if Y — ► K. is a vector bundle, there are the canonical splitting 
VY = Y xY and the morphism VY — > Y - so that s + s obeys the initial dynamic equation 
(]5aT) modulo the terms of order > 1 in s. 



3 The covariant Lyapunov tensor 

The collection of coefficients 

l k = d 3l k (6) 



of the variation equation (pb[) is called the Lyapunov matrix. Clearly, it is not a tensor 
under bundle coordinate transformations of the fibre bundle Y ([!]). Therefore, we introduce 
a covariant Lyapunov tensor as follows. 

Let a fibre bundle Y —>■ K be provided with a Riemannian fibre metric g, defined as a 
section of the symmetrized tensor product 

VV*Y^Y (7) 
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of the vertical cotangent bundle V*Y of Y — ► R. With respect to the holonomic coordinates 
(t,y k ,y k ) on V*Y, it takes the coordinate form 

a = ^9ij(t,y k W vdy*, 

where {cfo/} are the holonomic fibre bases for V*Y. 
Given a first order dynamic equation 7, let 

V*7 = l~d^y k W, ^' = Jr. (8) 

be the canonical lift of the vector field 7 (||D onto V*Y . It is a connection on V*Y — > M. 
Let us consider the Lie derivative L 7 g of the Riemannian fibre metric g along the vector 
field V*-y (|). It reads 

L i;j = (Lyfljy = + 7^5^ + (9 i 7 fe 5'fc i + dj~i k g ik . (9) 

This is a section of the fibre bundle (0) and, consequently, a tensor with respect to any 
bundle coordinate transformation of the fibre bundle ([]]). We agree to call it the covariant 
Lyapunov tensor. If g is an Euclidean metric, it comes to symmetrization 

Lij = <)r' + d,Y = U + 1/ 

of the Lyapunov matrix (|6]). 

Let us point the following two properties of the covariant Lyapunov tensor. 

(i) Written with respect to the initial date coordinates in Proposition [L], the covariant 
Lyapunov tensor is 



Jab 



d t g a b- 



(ii) Given a solution s of the dynamic equation 7 and a solution s of the variation 
equation (|5b|), we have 



The definition of the covariant Lyapunov tensor (^) depends on the choice of a Rieman- 
nian fibre metric on the fibre bundle Y. 

Proposition 2. If the vector field 7 is complete, there is a Riemannian fibre metric on Y 
such that the covariant Lyapunov tensor vanishes everywhere. 
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Proof. Let us choose the atlas of the initial date coordinates in Proposition Using 
the fibration ( : Y — > Y t=0 , one can provide Y with a time-independent Riemannian fibre 
metric 

g ab (t,y c ) = h(t)g ab (0,y c ) (10) 

where g a b{0, y c ) is a Riemannian metric on the fibre Y t= o and h(t) is a positive smooth 
function on R. The covariant Lyapunov tensor with respect to the metric ( [TOD is 

L a b = d t hg ab . 

Putting h{t) = 1, we obtain L = 0. □ 



4 The Lyapunov stability of a first order dynamic equation 

With the covariant Lyapunov tensor (||), we obtain the following criterion of the stability 
condition of Lyapunov. 

Recall that, given a Riemannian fibre metric g on a fibre bundle Y — * R, the instantwise 
distance p t (s, s') between two solutions s and s' of a dynamic equation 7 on F at an instant 
t is the distance between the points s(t) and s'(t) in the Riemannian space (Y t ,g(t)). 

Proposition 3. Let s be a solution of a first order dynamic equation 7. If there exists an 
open tubular neighbourhood U of the trajectory s where the covariant Lyapunov tensor ftQ) 
is negative- definite at all instants t > to, then there exists an open tubular neighbourhood 
U' of s such that 

lim [p t '{s,s') - p t (s,s')] < 

v— >oo 

for any t > to and any solution s' crossing U' . 

Proof. Since the condition and the statement of Proposition are coordinate-independent, 
let us choose the following chart of initial date coordinates in Proposition which cover 
the trajectory s. Put t = without a loss of generality. There is an open neighbourhood 
Uo C Yq n U of s(0) in the Riemannian manifold (Yq, g(0)) which can be provided with the 
normal coordinates (x a ) defined by the Riemannian metric g(0) in Yq and centralized at 
s(0). Let us consider the open tubular U' = C _1 (t^o) endowed with the coordinates (t,x a ). 
It is the desired chart of initial date coordinates. With respect to these coordinates, the 
solution s reads s a (t) = 0. Let 

s' a (t) = u a = const. 
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be another solution crossing U'. The instantwise distance p t (s, s'),t> 0, between solutions 
s and s' is the distance between the points (t, 0) and (t, u) in the Riemannian space (Y t , g(t)). 
This distance does not exceed the length 

1/2 



g ab (t,ru c )u a u b dT 



of the curve 

x a = TU a i r E [0, 1] 
in the Riemannian space (Y t ,g(t)), while 

p Q (s,s') =p (s,s') 



ill) 



(12) 



The temporal derivative of the function p t (s, s') ( pTT]) reads 

i 



dtp t (s,s') 



d t g a b(t,Tu c )u a u b dT. 



(13) 



Since the bilinear form d t g a b = L a b, t > 0, is negative-definite at all points of the curve 
(|T2|), the derivative (p~3|) at all points t > t Q is also negative. Hence, we obtain 



p t>0 (s,s') < p t>0 {s,s') < p {s,s') 



Po(s,s'). 



□ 



Corollary 4. The solution s obeying the condition of Proposition [| is Lyapunov stable 
with respect to the Riemannian fibre metric g. 

It is readily observed that Proposition [| states something more. One can think of 
the solution s in Proposition ^| as being isometrically Lyapunov stable. Of course, being 
Lyapunov stable with respect a Riemannian fibre metric g, a solution s need not be so 
with respect to another Riemannian fibre metric g', unless g' results from g by a time- 
independent transformation. 

Proposition 5. For any first order dynamic equation defined by a complete vector field 7 
(j^j on a fibre bundle Y — > R ; there exists a Riemannian fibre metric on Y such that each 
solution of 7 is Lyapunov stable. 

Proof. This property obviously holds with respect to the Riemannian fibre metric (|i0|) in 
Proposition ^| where h = 1. □ 
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Proposition |5| can be improved as follows. 

Proposition 6. Let X be a real number. Given a dynamic equation 7 defined by a complete 
vector field 7 there is a Riemannian fibre metric on Y such that the Lyapunov spectrum 
of any solution of 7 reduces to X. 

Proof. Recall that the (upper) Lyapunov exponent of a solution s' with respect to a 
solution s is defined as the limit 

K(s,s')=\\m±H Pt (s,s')). (14) 

t^co t 

Let us provide Y with the Riemannian fibre metric (|T0| ) in Proposition |2| where h = exp(At). 
A simple computation shows that the Laypunov exponent (|14]) with respect to this metric 
is exactly A. □ 

If the upper limit 

lim K(s, s') = X 

pt=o(s,s')—>0 

is negative, the solution s is said to be exponentially Lyapunov stable. If there exists at 
least one positive Lyapunov exponent, one speaks about chaos in a dynamical system ||. 
Proposition |6] shows that chaos in smooth dynamical systems can be characterized in full 
by time-dependent Riemannian metrics. 
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